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(54) Method for estimating or simulating parameters of a stratum structure 

(57) The present invention provides for a method for 
estimating or simulating one or more parameters of a 
stratum structure, comprising the steps of: 

determining a first model of the stratum structure; 
determining one or more variograms of one or more 
parameters in one or more of the strata; 
measuring the value of a parameter in a number of 
points; 

estimating on the basis of a variogram a probability 
density function (PDF) of the parameter at locations 
in the strata close to or at some distance from the 
locations where the parameter has been measured; 
and 

determining by using seismic data an acceptable 
value of the parameter at the location where the 
probability density function of the parameter has 
been determined using a variogram. 
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Description 

Particularly for locating oil and/or gas in the earth's 
crust it is important to obtain the best possible estima- 
tion of a stratum structure, for instance using seismic 5 
inversion techniques, whereafter it can be assessed 
whether (test) drilling for oil and/or gas can be carried 
out at a particular location on the earth's surface. Also 
after carrying out one or more test drillings, the costs of 
which are very high, examination of the stratum struc- 10 
ture continues, likewise using the measurement data 
obtained from the wells. 

In the article by Andr6 Haas et al: "Geostatical 
Inversion • A Sequential Method of Stochastic Reservoir 
Modelling constrained by Seismic Data", First Break vol is 
12, No 1 1. November 1994/561. geostatical simulations 
are used at an earlier stage of the inversion process. 
This know computer scheme is however quite compli- 
cated, e.g. because it is necessary to calculate syn- 
thetic traces. 20 

The present invention provides a method for esti- 
mating or simulating one or more parameters of a stra- 
tum structure, comprising the steps of: 

determining a first model of the stratum structure; 25 
determining one or more variograms of one or more 
parameters in one or more of the strata; 
- measuring the value of a parameter in a number of 
points; 

estimating on the basis of a variogram a probability 30 
density function (PDF) of the parameter at locations 
in the strata close to or at some distance from the 
locations where the parameter has been measured; 
and 

determining by using seismic data an acceptable 35 
value of the parameter at the location where the 
probability density function of the parameter has 
been determined using a variogram. 

According to the present invention geostatisticai 40 
techniques are combined with deterministic techniques 
which make use of seismic data, whereby the uncer- 
tainty in the solution of the model can be considerably 
reduced. 

The stratum structure is usually built up of sand- 45 
stones, shale and/or limestone-type rocks. Significant 
parameters are for instance the specific density, the 
porosity, the permeability and the like of such strata. 
The first model of the stratum structure on which itera- 
tive (better) estimations are performed is determined for so 
instance using seismic, geological and/or well data. The 
variogram of the parameter can be determined on the 
basis of the experience of a geologist, measurements 
on a well and/or from samples obtained from a well. 

Further advantages, features and details of the 55 
present invention will be elucidated with reference to the 
annexed drawings, in which: 



fig. 1 shows a diagram elucidating a preferred 
embodiment of the method according to the present 
invention; 

fig. 2 shows a graph of a variogram for use in the 
preferred embodiment according to the present 
invention; 

fig. 3 shows a graph of an example of a probability 
density function for use in the method according to 
the present invention; 

fig. 4 shows a flow diagram explaining the preferred 
embodiment of the method according to the present 
invention; 

fig. 5 shows a flow diagram explaining a further pre- 
ferred embodiment of the method of the present 
invention; and 

fig. 6 shows a flow diagram of a preferred embodi- 
ment wherein the methods according to fig. 4 and 5 
are combined. 

In a schematically designated area U (fig. 1) of a 
stratum structure, in which for instance two boreholes V 
and W are arranged, the parameter, for instance the 
porosity, must be determined of a large number of 
points - x, Xj, xj and Xk are shown - while the property of 
the structure at locations v and w at the respective bore- 
holes V and W are known with a certain accuracy, for 
instance because measurements are carried out with 
probes along the boreholes and/or samples are taken. 

In a preferred embodiment of the present invention 
a structural (geometric) model is first created on the 
basis of seismic measurement data at the earth's sur- 
face. Reflection surfaces which are considered impor- 
tant are formed to a structure of closed strata, for 
instance by the computer program EARTH MODEL 
commercially available from applicant, wherein a geolo- 
gist can make determined choices. The program 
EARTH MODEL defines the geological layers as well as 
the micro-structure in such a layer. It will be apparent 
that according to the present invention other methods 
may also be envisaged to obtain an initial simulated 
stratum structure. 

Variograms (fig. 2) (and/or histograms) are subse- 
quently defined for a particular parameter, such as the 
porosity in determined layers. In a variogram such as in 
fig. 2 the correlation or variance is determined between 
two values of a parameter given at different spatial loca- 
tions, i.e. as a function of the distance between two 
locations. A variogram between two different parame- 
ters, tor instance between porosity and permeability, is 
referred to as a cross-variogram. At a determined dis- 
tance the variogram becomes horizontal, above which 
distance the spatial correlation amounts to zero. The 
variance therefore has a constant value at a distance 
beyond this distance. In the case a variogram is aniso- 
tropic in a particular stratum, this distance has a differ- 
ent value in each case for different directions, thus 
forming as it were an ellipsoid. 

In order to determine the parameter in the point x 1( 
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which can be chosen randomly (or according to more 
advanced probability techniques which are not dis- 
cussed here), a number of points, for instance a, b, c 
and d. where the value of this parameter is known (or is 
already present in the model and has been determined 
previously in the same manner) are chosen in the vicin- 
ity of this point x v Via the probability density function 
found on the basis of the variograms a value for the 
parameter is chosen randomly and stored for the point 
x v The thus simulated point is added to the known 
value and the procedure is repeated until all points xj in 
the space U have been simulated and a value assigned 
thereto. 

The above stated preferred method differs from the 
so-called Kriging method (or in the case of cross-vario- 
grams co-Kriging) in that therein the most probable 
value is determined from the variograms (see C. V. Deut- 
sch and A.G. Journel (1992), GSLIB: geostatistical soft- 
ware library and user's guide. Oxford University Press, 
New York). 

In the stochastic techniques an estimation of the 
uncertainty in the estimated model reservoir is availa- 
ble, wherein the reservoir models fit available geostatis- 
tical information, such as variograms and histograms 
(i.e. wherein the data has a Gaussian distribution). 

According to the present invention the models, such 
as obtained with the above mentioned stochastic simu- 
lation, are subsequently checked against seismic data. 
Instead of the parameters being estimated exclusively 
according to stochastic simulations, after obtaining a 
first stochastic model it is determined using seismic 
data whether or not a new iteration point fits better in the 
model than the previously found point. 

This takes place by convolving synthetic seismic 
data with a "wavelet" and comparing the outcome with 
the measured seismic data. If the value of the iteration 
point is no better than the previous point, this iteration 
step is rejected and the old value of the point is stored in 
the model. 

Due to the combination of stochastic simulation and 
application of seismic data, wherein iterations can con- 
tinue until a satisfactory correlation between the reser- 
voir model and the seismic data is obtained, there 
results an estimation or model for the values of the 
parameters in the various points which suffice optimally 
on the basis of the geostatistical data. 

In the flow diagram according to fig. 4 various 
aspects are elucidated further: firstly (10) a point or 
node Xf is selected randomly. The conditional local prob- 
ability density function of the reservoir property is then 
estimated (11); the realization is drawn therefrom (12), 
which procedure is repeated until all grid nodes are 
dealt with. 

In fig. 5 a similar method is shown wherein firstly a 
grid node (20) is selected at random, the conditional 
local probability density function is then determined 
subject to a particular lithology (indicator) (21 ). whereaf- 
ter the realization lithology from the local probability 
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density function is drawn (22) until the random path (23) 
is complete. 

Fig. 6 shows the preferred embodiment, wherein 
the procedures of fig. 4 and 5 are combined and 
wherein the convergence step for the seismic data is 
shown in block 31. 

First experiments with the method according to the 
present invention have shown that the cumulative prob- 
ability in the 10/90% probability interval is considerably 
smaller in the method according to the present invention 
than using the known geostatistical methods. 

The present invention is not limited to the above 
described preferred embodiment; the requested rights 
are determined by the following claims, within the scope 
of which many modifications can be envisaged. 

Claims 

1. A method for estimating or simulating one or more 
parameters of a stratum structure, comprising the 
following steps: 

- determining a first model of the stratum struc- 
ture; 

determining one or more variograms of one or 
more parameters in one or more of the strata; 
measuring the value of a parameter in a 
number of points; 

estimating on the basis of a variogram a proba- 
bility density function (PDF) of the parameter at 
locations in the strata close to or at some dis- 
tance from the locations where the parameter 
has been measured ; and 
determining by using seismic data an accepta- 
ble value of the parameter at the location where 
the probability density function of the parame- 
ter has been determined using a variogram. 

2. A method for estimating or simulating one or more 
parameters of a stratum structure, comprising the 
steps of: 

a) obtaining seismic data of the stratum struc- 
ture; 

b) determining from the obtained seismic data 
a first model of the stratum structure; 

c) determining from the first model one or more 
variograms of one or more parameters in one 
or more strata of the stratum structure; 

d) measuring a value for the one or more 
parameters at a plurality of points in the stra- 
tum structure; 

e) estimating on the basts of the one or more 
variograms a probability density function (PDF) 
of the one or more parameters at other points 
in the strata one of adjacent and distant from 
the points where the one or more parameters 
were measured; and 
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f) determining by using the probability density 
function a value for each of the one or more 
parameters in the strata at the points where the 
probability density function of the parameter 
were determined; s 

g) combining the values determined in step f) 
wtth the corresponding one or more parame- 
ters determined from the first model; 

h) obtaining from one or more parameters 
determined in step g) acceptable values for the 
one or more parameters in the stratum struc- 
ture; and 

0 generating from the acceptable values for the 
one or more parameters a second model of the 
stratum structure. 

3. A method as claimed in daim 1 or 2, wherein for 
each location two parameters are simulated using a 
cross-variogram. 

4. A method as claimed in daim 1 . 2 or 3, wherein the 
probability density function is selected on the basis 
of an indicator for the type of lithology to be 
expected. 

5. A method as claimed in claims 1 -4, wherein the one 
or more variograms are determined on the basis of 
available data of a stratum obtained from samples 
thereof. 

6. A method as claimed in daims 1 -5, wherein one or 
more parameters are measured in a borehole. 

7. A method as claimed in claims 1-6. wherein the 
steps are repeated until an estimation is obtained of 
desired accuracy. 

8. A method as claimed in claims 1 -7, wherein the one 
or more variograms are determined on the basis of 
available data of the stratum obtained from samples 
thereof. 

9. Trie method as claimed in claim 1 , wherein step (h) 
comprises the steps of: 



10. The method as daimed in claim 9, further induding 
repeating steps (g) - CO until a desired correlation is 
obtained between the one or more parameters in 
the first model and the one or more parameters in 
. the second model. 
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convoluting the one or more parameters deter- 
mined in step (g) with a wavelet to obtain con- 
voluted one or more parameters; 
comparing the convoluted one or more param- 
eters with corresponding one or more parame- so 
ters determined from the first model; and 
retaining the obtaining convoluted one or more 
parameters as the acceptable values if the 
obtained convoluted one or more parameters 
improve a cumulative probability of the proba- ss 
bility density function over the one or more 
parameters from the first model. 
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